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Abstract 

We calculate event rates and demonstrate the observational feasibility of very high energy 
muons (1 TeV-1000 TeV) in a large mass underground detector operating as a pair-meter. 

This energy range corresponds to surface muon energies of ~(5 TeV - 5000 TeV) and primary 
cosmic ray energies of ~ (50 TeV - 5 xlO 4 TeV). Such measurements would significantly assist 
in an improved understanding of the prompt contribution to v e ,v^ and g fluxes in present 
and future ultra-high energy neutrino detectors. In addition, they would shed light on the 
origin of and possible compositional changes at and around the observed ’knee’ in the cosmic 
ray spectrum. 

1 The Cosmic Ray Spectrum and the Knee 

Cosmic ray studies, with the spectrum extending over ten decades in energy, have proved to be 
fertile terrain for furthering our knowledge of both astrophysics and particle physics ( reviews may 
be found in [D121EJ H). They have traditionally provided us with clues for the existence of new 
particles and the physics associated with them, which have later been confirmed by detailed accel¬ 
erator experiments. In fact, prior to 1950 and the advent of modern accelerator technology, they 
provided the only means of studying high-energy particle production and interactions. Addition¬ 
ally, as a result of our attempts to undestand the origin of cosmic rays, they have contributed to 
our knowledge of acceleration via shocks, and the propagation of charged particles in the galaxy 
and heliosphere. 

The cosmic ray spectrum, characterised by a steeply falling power-law behaviour over its entire 
range, exhibits two transition regions where the slope changes noticeably: 

• A steepening of the spectrum occurs around E fa 5 x 10 6 GeV, i.e. the index 7 describing 
the power-law behaviour of the differential flux, dN/dE ~ E 1 , changes from 7 m —2.7 to 
7 as —3.1; leading to the feature called the ‘knee’. 
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• A flattening of the spectrum occurs around E « 5 x 10 9 GeV, i.e at the “ankle”; with the 
index 7 changing back to ~ 2.4 — 2.7. Beyond the ankle, in the realm of ultra high energy 
cosmic rays, data u u u\ is sparse and conflicting, but highly intriguing. While we will 
not address the interesting puzzle in this regime here, a discussion of the various issues may 
be found in [HI El, and a recent assesment of the shape of the spectrum based on current 
knowledge can be found in m- 

The physical reason for the existence of the knee is at present an unresolved problem of great 
significance to understanding the origin of galactic cosmic rays. It is generally believed that the 
reasons underlying this distinctive shift in the spectrum are astrophysical in nature, as opposed 
to those stemming from a change in hadronic interactions at these energies which, at present, are 
not within the reach of existing accelerators. This conclusion is based on the observed correspon¬ 
dence between independant measurements of the muon number spectrum, Cerenkov radiation and 
hadronic constituents of air-showers mmrnmm- While the reasons for the shift in the 
spectrum remain un-understood, these data exhibit an expected co-relation which supports the 
absence of radically different physics interactions at these energies. 

While the case for the existence of new physics being at least partially responsible for a shift in 
the spectral index is not wholly without motivation 3 * , we stress that it appears unlikely that this can 
be empirically corraborated or refuted in the near future in CR measurements. This is because, as 
we shall discuss below, uncertainties (in the knee region) in the CR compositon and prompt muon 
and neutrino contributions would likely overshadow evidence of such new interactions. In any case, 
this hypothesis will be thoroughly probed in the near future by the Large Hadron Collider (LHC) 
at CERN which will operate at a center of mass energy of ^/s = 14 TeV. 

2 Uncertainties in the Muon and Neutrino Fluxes in the 
Knee region and beyond 

As stated earlier, present data [20] when culled and correlated, appear to favour one or more 
astrophysical reasons for the existence of the knee. These include it being a rigidity-dependant 
effect (originally proposed in jT[]) related to the (different) maximum acceleration energies for 
different nuclei either in the cosmic ray source itself or during the propagation process. Data from 
surface air-showers and optical detectors indicate, without being conclusive, that the average mass 
of the cosmic ray spectrum nuclei differs before and after the steepening at the knee. In particular, 
there appears to be some evidence [221125 j that the composition is heavier above the knee region. If 
this is true, then, as discussed in 1231. significant suppression of the very high energy (> 10 5 GeV) 
muon and neutrino fluxes resulting from CR interactions in the atmosphere and in the interstellar 
medium can occur. 

A second major factor in determining the enhancement (or lack thereof) of muon and neutrino 
fluxes above several TeV are the uncertain magnitudes of the prompt (i.e. those resulting from 
heavy meson decay, notably, charm mesons and heavier composites) fluxes of both. At low (i.e. 

3 This region in E p , (i.e, the energy of the primary particle) corresponds to several TeV in center of mass 
energies. Thus there are many conjectures for physics beyond the Standard Model which come into play, e.g. 

SUSY, technicolour, large extra dimensions etc. These could lead (via new particle production and decay) to energy 
being channelled into muons, neutrinos or other secondary particles in a manner that present cosmic ray experiments 

are insensitive to, causing the shift in the (measured) energy spectrum dm HH HH HD 
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~ GeV) energies, the cosmic ray induced neutrino and muon fluxes receive their dominant contri¬ 
butions from the decays of 7r and K mesons, whose interaction lengths significantly exceed their 
decay lengths m ESI EZ1 EE] • ( These fluxes are henceforth referred to as the conventional fluxes 
in what follows.) This sitiuation changes at ~ TeV energies, and secondary interactions of these 
particles become possible, leading to the production of heavy short lived hadrons. While upper 
bounds on the flux of muons and neutrinos have been provided by several experiments e.g LVD 
m, AKENO m an d AMANDA jO], they still allow for a very large possible range of prompt 
flux magnitudes. 

Present phenomenological predictions for the diffuse fluxes of these prompt muons and neutrinos 
can differ by about two orders of magnitude E21 53 E3I ESI EZj The sources of this large 
uncertainty lie, to a significant extent, in the choice of charm production models. For instance, 
differing predictions arise from models based on perturbative QCD (pQCD) with a K factor [ IK] , 
next-to-leading order (NLO) pQCD j-bfit I54j , quark-gluon string models and recombination quark- 
parton modelsP2J etc. In general, QCD based models must contend with a large uncertainty 
associated with the extrapolation of the gluon parton distribution function g(x) to small fractional 
momentum x < 10 -5 . Theoretical models generally assume 

xg(x) ~ x~ x , 

where A is in the range 0 — 0.5, and fluxes depend strongly on the chosen value of A. We note 
that depending on the model, the prompt muon and neutrino fluxes from charm decay exceed 
the corresponding conventional fluxes (from 7 r and K decays) somewhere between (surface) muon 
energies of few tens of TeV and few PeV [SB]. Reliable measurements of muon fluxes in this range 
would thus, at the very least, help in establishing the reliability of a particular class of models. 

While we do not give a detailed account of the flux predictions from all the different models, 
we attempt to give a representative idea in our calculations of the variation possible even within 
a given charm production model. While the conventional muon flux from 7 r and K decays is well 
understood and fairly firm, the prompt flux predictions are subject to variations resulting from 
different parton distribution functions and choices of the factorization and renormalisation scales 
as mentioned above. We thus stress the need for better empirical determination of the muon (and 
associated neutrino fluxes) in this region, a topic which is elaborated upon in the next section. 

3 The significance of measurements of the muon and neu¬ 
trino fluxes in the knee region 

With very few exceptions, available data on muons above several TeV comprise of measurements 
of the number spectrum rather than the energy. This is primarily due to the size and density 
requirements imposed on detectors by the significant penetration lengths acheived by high energy 
muons. 

The desirability of improved and statistically significant muon energy measurements in the few 
TeV to few hundred TeV region stems from (at least) three reasons: 

• As mentioned above, the physical origin and composition of cosmic rays in this energy range 
is currently obscured by a paucity of data on VHE muon and neutrino fluxes. Observations 
would thus certainly illuminate the current debate on the reason for the occurance of and 
compositional changes at the knee. 
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• Also, as discussed above, QCD related theoretical uncertainties dominate the predictions 
of the prompt contribution to muon and neutrino fluxes. As emphasized in [ 551155] , the 
measurement of down-going muon fluxes would provide a valuable handle in their reduction. 

• Both the conventional and prompt muon fluxes at these energies are closely related to the 
associated neutrino fluxes. For prompt contributions, this is is because the kinematics of 
charmed particle decay and the corresponding semi-leptonic branching ratios ensure that the 
v e and Vp fluxes are identical upto a few percent to the muon fluxes in this energy range, 
regardless of the choice of the charm production model or of A. The conventional neutrino 
flux, on the other hand, is about 10 % of the conventional muon flux. At the energies of 
interest here, neutrinos resulting from cosmic ray interactions in the atmosphere and in the 
inter-stellar medium constitute the most important background to searches for diffuse fluxes 
of ultra high energy (UHE) neutrinos [321 3 dH S3 ES. S3 EE] from cosmological sources 
(e.g. active galactic nuclei, gamma-ray bursts etc.) in neutrino telescopes like AMANDA 
m, ICECUBE m\ an d NEMO [Hj. Thus, they are an important obstacle to the much- 
anticipated detection of such energetic point sources in these detectors. Empirical data on 
downgoing muons from cosmic rays would prove invaluable in understanding this background 
to the UHE neutrino signal, since it can be co-related to the neutrino flux. 

In the context of the points above, it is relevant to stress the importance of being able to dis¬ 
entangle the prompt (due to the decay of produced heavy mesons) and conventional (from n and 
K decays) and diffuse UHE (from extra-galactic sources) contributions to the neutrino fluxes. 
Methods to enable this have been studied, and are based on the differences in zenith angle, depth 
and spectral dependance between these fluxes |421 53] and on the isolating capability of showering 
u e charged current events via a break in the spectrum [ 55] , 

Having emphasized the importance of muon energy measurements in the several TeV to several 
hundred TeV range, we proceed in the next section to study the potential of the pair meter 
method 123123123 as applied to such measurements made in a large iron calorimeter (50 kT) 4 . 
Since individual muon energies will become measurable using this technique, it will be possible 
to augment the sparse existing data on cosmic ray muons in the important range where they 
have surface energies of 5 — 5000 TeV. Furthermore, these observations can be combined with 
balloon-based experiments ( e.g TRACER [2Q]) and upcoming hybrid air-shower experiments (e.g 
KASCADE-Grande J20J and LOPES pH] ) to enhance our understanding of the issues discussed 
above. We mention here that this range in muon surface energy roughly corresponds to a range of 
50 — 5 x 10 4 TeV in primary cosmic ray energy, which is crucial to an enhanced understanding of 
the origin of the knee. 

In the remainder of the paper, we first provide a discussion of the pair-meter technique and 
the pair production cross section which results in the observed cascades. This is followed by 
a brief description of a typical large-mass iron calorimeter. Subsequent to this we summarize 
the interactions and losses of muons in matter enroute to an underground detector, and their 
incorporation into our calculation. We then calculate anticipated event rates for a 50 kT detector 
and demonstrate that even after accounting for energy losses in the surrounding rock, event rates 
can be appreciably large for the 1 — 1000 TeV range, corresponding to surface muon energies in 
the range of several TeV to several PeV. 

4 Such a detector is currently being planned for location in India [281 
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Figure 1: Differential cross section vda/dv vs. v ^inverse of the relative energy transfer)for pair 
produetionisolidl jHH]. bremstrahlungfdotted) [HTj and photonuclear (dashed) [H2] processes. 

4 The Pair Meter method and the associated Pair Produc¬ 
tion Cross Section 

Due to the penetrating power of muons, their energy measurements require techniques which differ 
from those employed for photons, hadrons and electrons. Furthermore, muon energy measurement 
methods which work well in the GeV range (magnetic spectrometry or measuring Cerenkov radia¬ 
tion) are rendered impractical in the TeV range primarily due to requirements of size imposed by 
the combination of high energies and a steeply falling spectrum. 

The pair meter technique [23 123 I2Z| skirts some of the disadvantages of traditional muon 
detectors by relying on a somewhat indirect method, i.e. the measurements of the energy and 
frequency of electron-positron pair cascades produced by the passage of a high energy muon in 
dense matter. A reliable reconstruction of the muon energy in this method is based on the following: 

• The cross section for e + e” pair production by a muon with energy E with energy transfer 
above a threshold E 0 grows as ln 2 (2m e E^/m ll E 0 ), where and m e are the muon and 
electron masses respectively. 

• Defining v = Eq/'E^, above v = 10, this cross section dominates those for other muon energy 
loss processes which generate observable cascades in its passage through dense matter, e.g. 
/i — N inelastic scattering and bremstrahlung emission. This is demonstrated in Figure ^ 
where we compare the differential cross sections for these various interactions as a function 
of v~ l . 

• The energy lost to each cascade resulting from e + e~ pair production is a very small fraction 
(about 1(T 2 ) of the muon energy for the range of u -1 which we focus on here. 

• The dependance of the pair production cross section on E^/E 0 then allows one to infer the 
muon energy by counting the number of interaction cascades N in the detector with energies 
above a threshold E 0 . 
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We now make the above statements more precise. In the approximations 



2 m e 


m 




and 

E 0 > 2m e 189 v / e^" 1/3 ~ 0.3Z- 1/3 GeV, 


(Z—atomic nnmber= 26, for iron) both of which are valid for the choice of E 0l E ]t for which we 
present results below, the expression for the differential pair production cross section is given bv|59| 


do 14a (Km e En\ 

"a," 9 ■ 

where a = 1/137 and k ~ 1.8. to is the radiation lenth (r.l) which is given by 


to = ( — + — olt^Na ln(189Z 1 / 3 ) > ] 

V ) 


( 1 ) 

( 2 ) 


Here Aw is the gram atomic weight (for iron, this is 56 grams), is the classical electron radius 
and Na the Avogadro number. For iron, this gives to = 13.75 gm/cm 2 . 

The average number of interaction cascades M above a threshold Eq for v < 10~ 3 is given by 


M(E 0 ,E,) — Tt 0 o(E 0 , Efj) , 


(3) 


where T is the thickness of the target in units of t 0 and o(E 0 , E tI ) is the integrated cross section 
(in units of cnr/grn), 


o(E 0 , E^) ~ 


7a 


In 


2 ( KrrieEf, 


97rt 0 \ V E 0 m^ 


+ C 


(4) 


where C ~ 1.4. 

At this point it is important to mention that the capability and effectiveness of the pair meter 
method for high energy muons has been tested and demonstrated by the NuTEV/CCFR collabora¬ 
tion, as described in [HO] • The calculations which follow are performed for a 50 kT iron calorimeter. 
Our prototype is based on the suggested design for INO; see J2B, f° r details. The dimensions of a 
50 kT detector of this type would correspond to (approx) 15 m x 15 m x 45 m. A muon traversing 
a 20 m path in this detector corresponds to a path-length of ~ 1145 r.l. In what follows, we assume 
a (conservative) “average” path-length of 1000 r.l for the typical muon and calulate the number of 
observable cascades produced by it, for different cascade thresholds and muon energies. Figure El 
shows the average number of cascades above a threshold energy Eq produced by a muon entering 
the detector with energy E^ and T = 1000 r.l.; for three different choices of Eq, i.e. 1 GeV, 10 
GeV and 100 GeV. Quantitatively, we note that this leads to a E^ = 100 TeV muon generating 
approximately 40 cascades, each of energy greater than E 0 = 10 GeV and 10 cascades with energy 
in excess of 100 GeV. By counting the cascades for several choices of thresholds for a traversing 
muon, one obtains a reliable estimate of its energy. 

It is also relevant to remark here that the relative energy measurement error, 5E^/E^ in the 
pair meter is given by 


SEp/En = 


' 97T 
28 aT 



(5) 
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Muon Energy (GeV) 

Figure 2: Average number of cascades above a threshold E 0 vs. muon energy for E 0 — 1 GeV 
(solid line), 10 GeV (dotted) and 100 GeV (dashed), with T fixed to 1000 r.l. 

For v = (10 -3 — 1CT 2 ), which is the range we focus on here, this allows a liberal tolerance for error 
in the measurements of individual cacade energies. We note also that the errors do not worsen 
with increasing muon energy, which is an important advantage of the pair-meter technique. 


5 The Surface Muon Energy determination for Underground 
Events 

It is important to co-relate the measured muon energies in an underground detector to their surface 
energies, which we take to be those that would be observed were our detector placed on the surface 
of the earth. This requires a calulation of the energy loss as the muon traverses the rock between 
the earth’s surface and the detector. 

These losses originate from ionization, bremsstrahlung, pair production and photonuclear in¬ 
teractions. They can be effectively parametrized oTi lUT , 52; for E t , > 1 TeV, since the average loss 
increases predominantly linearly with energy, 

= < 6 > 

where a parametrizes the contribution from ionization of muons and f3 encapsulates the contribu¬ 
tion from bremsstrahlung, pair production and photonuclear processes. Note that a and /3 carry a 
very weak(intrinsic) energy dependance. It is thus appropriate to assume that the average muon 
energy at depth X is 

(E,XX)) = (b; + ^ e-^ , (7) 

where E^ is the initial surface muon energy. One may use this to write down the minimum surface 
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energy required of a muon to reach a depth A" as, 


^min 



( 8 ) 


From Eq. 0 we get the relation between initial energy E s and degraded energy of muon E /t 
after travelling a distance X as, 


E 


s 



JX 


a 

1 ' 


The differential muon flux at a depth X is given by, 

dN _ dN 
dEdE s ^ 


where is the initial muon flux with surface muon energy E s f] . 


(9) 


( 10 ) 



Figure 3: Degraded muon energy, E vs. surface muon energy, E s after passing a depth of 3.5 x 
10 5 g/cm 2 , corresponding to the proposed INO location. 

The a and f3 in the analytical expression can be obtained for standard rock. The depth relevant 
to the INO detector’s proposed location is 3.5 x 10 5 gm/cm 2 . The values are, 

/3 = 4 x 10 -6 cm 2 /g, = 675 GeV. (11) 

P 

Figure 0 shows the degraded muon energies (he. those measured for muons entering the detector 
after traversing the rock) vs their corresponding surface energies after losses are accounted for in 
the manner described above. We note that typically, E^ ~ (2 — 5) x E^. 





Flux model 

PDF 

Scales 

a 

b 

c 

d 

PRS1 

CTEQ3 

M = p, = m c 

5.37 

0.0191 

0.156 

0.0153 

PRS 2 

CTEQ3 

M = 2/i = 2 m c 

5.79 

0.345 

0.105 

0.0127 

PRS3 

D 

M — 2jl — 2 m c 

5.91 

0.290 

0.143 

0.0147 


Table 1: PRS parameters for the prompt muon and antimuon fluxes. m c is the mass of the charm 
quark. 



Number of muons per solid angle entering the detector in 5 years 

E^TeV) 

conv+TIG 

conv 

TIG 

PRS1 

PRS2 

PRS3 

1 

1.035 x 10 7 

1.03 x 10 7 

37461 

55482 

95489 

136871 

10 

52486 

51282 

1204 

2952 

5341 

10443 

50 

770 

696 

74 

236 

431 

1104 

100 

127 

106 

21 

73 

134 

387 

200 

22 

16 

6 

22 

40 

129 

300 

8 

5 

3 

11 

19 

66 

400 

4 

2 

2 

6 

11 

41 

500 

2 

1 

1 

4 

7 

28 

600 

1.5 

1 

.5 

3 

5 

20 

700 

1 

.5 

.5 

4 

7.5 

31 

800 

.8 

.35 

.5 

1.5 

3 

12 

900 

.65 

.25 

.37 

1.25 

2.5 

10 

1000 

.5 

.2 

.3 

1 

2 

4 

10000 

.0025 

.0003 

.0022 

.007 

.013 

.08 


Table 2: Number of muons per solid angle entering the detector over 5 years for various energies 
of the entering muon, E [L . 


6 Muon Fluxes 


Extensive predictions and studies lifil liTl li£l 1521IH3115il I5H1 IKfil 157] for prompt cosmic ray muon 

fluxes at very high energies exist in the literature, as mentioned earlier. For our representative 
calculations of muon event rates, we have used the relatively conservative predictions for charm 
induced fluxes given in HUES]. The large variation in muon rates possible due to flux uncertainties 
even when these fluxes are used is amply reflected in our results, most noticeably in Table 2. One 
would expect much larger variations if the full range of prompt flux models available is used to 
calculate event rates. 

In [48] (henceforth referred to as the TIG flux), the conventional and prompt fluxes have been 
parametrized as 

dN N 0 E -T '- 1 , . 

-= — - (12) 

dE l+AE { ’ 


for E < E a . and as 

dN N'qE-^'- 1 
dE ~ l + AE 


(13) 


for E > E a . For the conventional muon flux N 0 = 0.2, N' 0 = 0.2, 7 = 1.74, 7 " = 2.1, E a = 5.3 x 10 5 , 
A = 0.007. 
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Number of cascades per muon for different thresholds Eq in GeV 


E, 

e; 

5 

10 

20 

50 

100 

300 

500 

1000 

5000 

1 

6.1 

3.08 

2.56 

3.78 







10 

40.26 

17.28 

10.99 

6.43 

3.08 

2.56 





20 

83.16 

25.3 

17.28 

10.99 

5.34 

3.08 





50 

205 

38.58 

28.26 

19.67 

10.99 

6.43 

2.78 

2.56 



100 

407.58 

50.63 

38.58 

28.26 

17.28 

10.99 

4.58 

3.08 

2.56 


200 

813 

64.43 

50.63 

38.58 

25.30 

17.28 

8.11 

5.34 

3.08 


300 

1218 

73.3 

58.49 

45.42 

30.8 

21.76 

10.99 

7.46 

4.19 


400 

1624 

79.96 

64.43 

50.63 

35.06 

25.3 

13.39 

9.33 

5.34 


500 

2029 

85.33 

69.24 

54.89 

38.58 

28.26 

15.45 

10.99 

6.43 

2.56 

600 

2435 

89.85 

73.3 

58.49 

41.58 

30.8 

17.28 

12.47 

7.46 

2.58 

700 

2841 

93.76 

76.83 

61.64 

44.21 

33.05 

18.91 

13.82 

8.43 

2.6 

800 

3246 

97.23 

79.96 

64.43 

46.56 

35.06 

20.4 

15.06 

9.33 

2.72 

900 

3652 

100.33 

82.77 

66.95 

48.69 

36.9 

21.76 

16.21 

10.18 

2.89 

1000 

4057 

103.16 

85.33 

69.24 

50.63 

38.58 

23.02 

17.28 

10.99 

3.08 

10000 

40554 

174.84 

151.24 

129.38 

103.16 

85.33 

60.63 

50.63 

38.58 

17.28 


Table 3: Number of cacades above thresholds E 0 = 5,10, 20, 50,100, 300, 500,1000, 5000 GeV per 
muon. Here E ' M is the energy of the muon in TeV entering the detector, and is its corresponding 
energy in TeV at the surface of the earth, assuming it traversed a depth of rock corresponding to 
3.5 x 10 5 gm/crn 2 . 


For the prompt muon flux N 0 = 1.4 x 10 5 , N' 0 = 4.3 x 10 4 , 7 = 1.77, y' = 2.01, E a = 9.2 x 10 5 , 
A = 2.8 x HT 8 . 

The second set of representative prompt muon fluxes we use are calculated in [51] (henceforth 
referred to as the PRS1,PRS2 and PRS3 fluxes). The differences in the three fluxes originate in 
different choices of parton distribution functions(PDF) and factorisation (M) and renormalisation 
scales(/i) of the theory. These fluxes can be convieniently parametrized |3! as follows 

dN _ -t r^~ a -\-b X + C x 2 —d X 3 (~\A\ 

dE~ ’ K } 

where x = Logio(E/ GeV), with a, b, c and d as in Tabled 

In Figure 0 we show the conventional (TIG) and prompt (TIG and PRS) surface muon fluxes. 
Uncertainties in the conventional flux, unlike the prompt case, are not major, hence we have shown 
only the TIG parametrization. We note that depending on the flux model, the prompt fluxes rise 
above the conventional flux for (surface) muon energies between 200 TeV and 1000 TeV. In terms of 
(degraded) muons entering the detector, we see from Figured and Figured that this corresponds 
to measured muon energies of several tens of TeV and several hundreds of TeV. Thus, we note 
that underground muon measurements in this range will help reduce the present uncertainties 
in deducing the charm contributions to muon and neutrino fluxes. Our calculations provide a 
quantitative estimate of the feasibility and the potential of these measurements to accomplish this. 
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E^GeV) 


Figure 4: (E^ 3 x Flux) vs. surface muon enegy E* for total TIG flux(solid), conventional (dotted), 
TIG prompt (short dashed), and prompt fluxes corresponding to PRSl(short spaced dots), 
PRS2(large dashed) and PRS3(large spaced dots). 



E u (GeV) 


Figure 5: (E 3 x Flux) vs. energy of muon entering the underground detector E M , ( for the flux 
models listed in the previous figure) after passing through the rock distance of 3.5 x 10 3 gm/cm 2 . 

. Labelling of plots is identical to figure 
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E^GeV) 

Figure 6: Number of muons entering the 50 kT detector in 5 years per solid angle vs. muon energy. 
Energy losses in the surrounding rock are taken into account. Depth of the detector is assumed to 
be 3.5 x 10 5 g/cm 2 . 

7 Results and Discussion 

We are now in a position to calculate the expected cascade events for a 50 kT detector in the 
energy range of interest discussed above. While an entering muon in this energy range will produce 
observable cascades, the number entering the detector over a given period is limited by the sharply 
falling fluxes at these energies. It is thus pertinent to obtain a quantitative measure of this by 
estimating n lt , the number of muons above a given threshold entering the detector per ster-radian 
for an exposure of t years, 

r °o gL/V 

= / dE fl ——A x t , (15) 

J Eth 

where A is the exposed area of a 50 kT iron detector. This is shown in Figure |H1 and Table [21 
We note that while the number of entering muons for the lowest energy in Table [21 k e. 1 TeV 
is very large, one also obtains an observable number, i.e. 1-3 events after integrating over solid 
angle (considering that there is no “back-ground” as such for such events) over the 5 year period 
even for E = 1000 TeV for the most conservative flux choice (TIG). These energies delineate the 
muon energy range accesible. The number of entering muons for all choices of PRS fluxes will be 
substantially higher, as shown. Even for the most conservative (TIG) flux choice, one expects good 
observational capability upto several hundred TeV 5 . 

The number of cascades per muon above E 0 = 5, 10, 20, 50, 100, 300, 500, 1000, 5000 GeV 
respectively using Eq. [21 and El are tabulated in Table [21 Table [21 also lists the surface muon energy 
E* corresponding to the underground muon energies E fl . While these are sample choices, it is clear 

5 We have given the per-steradian rates here. In order to predict the rate integrated over angle, predictions 
must carefully take into account the depth dependance of the fluxes for a particular detector and its surrounding 
topography. 
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that they can be further optimized based on the muon energy that one wants to observe to good 
statistical accuracy. 

The total number of cascade events per ster-radian N C (E 0 ) (above a given threshold E 0 and for 
a 50 kT x 5 yr exposure) is given by 

N C (E 0 ) = M (#„, Eq)^ , (16) 

where M(E fl , E 0 ) is the cascade number calculated above in Section 2.1. From Tableland Table 
01 we observe that this number is considerable for most thresholds, promising rich observational 
capabilities. For example, for the (most conservative) conventional+TIG flux model which we use 
as our benchmark, one finds that at even at muon energies of 1000 TeV,one can produce 51 events 
per solid angle for a threshold of 5 GeV and ~1 event per solid angle for a threshold of 5000 GeV. 
Expectedly, the PRS models predict significantly more events compared to these estimates. 

We have used the TIG flux as a benchmark to establish observability, since it leads to the most 
conservative event rate predictions. All other flux parametrisations lead to higher predictions. We 
note that even though TIG and PRS are not vastly different from each other in a qualitative sense 
since both are based on perturbative QCD inputs, their event rates in a large mass pair meter differ 
significantly. Indeed, the variations amongst fluxes in the same family(PRSl, PRS2, PRS3) are also 
large. Thus, the muon event rate can act as a soft ( i.e not definitive, given the large uncertainties 
in the QCD predictions) discriminator between various prompt flux models and provide pointers 
to the physics input that should guide their development. Similarly, this rate provides a tool to 
better understand the present spectral uncertainities in the cosmic ray knee origin. 

8 Conclusions 

Our main results are presented in Figs H3 and 01 and Tables El 01 From these we (conservatively) 
conclude that underground muon energy measurements for an energy range of E M of 1-1000 TeV 
are possible with a 50 kT iron detector 6 running for 5 years. This will enable a better handle on 
the very high energy muon fluxes between several TeV to about 5 PeV, and consequently illuminate 
our estimates of the background muon and neutrino fluxes for ultra high energy neutrino detectors 
and lessen present uncertainties in charm production models. As emphasized earlier, the prompt 
muon flux is a measure of the prompt v e and flux, hence its importance to ultra high energy 
neutrino astronomy cannot be underestimated. 

The observable muon energy range discussed in our results also corresponds to a range of 50 
TeV to 50 PeV in primary cosmic ray energies. This range is crucial to an understanding of the 
origin of knee and our calculations demonstrate the feasibility and potential resulting from muon 
measurements for a better understanding of the origin of the knee. 

A detailed and comprehensive set of predictions for a given large mass detector necessarily 
requires a much more elaborate calculation of the muon losses than what is presented here, since 
local topography plays an important role in determining the surface muon energy corresponding 
to a measured muon energy. Our aim in this paper has been more to demonstrate observational 
feasibility rather than to make precise predictions. Thus the calculations here show that very high 
energy muon measurements are possible in a large iron calorimeter and can aid in illuminating 

6 As stated earlier, we have used a depth corresponding to the proposed INO site for specificity, however, the 
results can be easily generalised to other depths, as should be obvious. 
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three important outstanding questions which address partially overlapping isuues, one in cosmic ray 
physics, the second in theoretical QCD and and the third in ultra high energy neutrino astronomy. 

Acknowledgment: We thank Namit Mahajan for helpful discussions. RG would also like to 
thank Francis Halzen for informative discussions and John Beacom for pointing out some helpful 
references. 
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